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ABSTRACT 



Evaporative evolution of stellar clusters is shown to be relaxation limited when 
the number of stars satisfies N » N^, where — 1600. For a Maxwell velocity 
distribution that extends beyond the escape velocity, this process is bright in that 
the Kelvin-Helmholtz time scale, f^^treiax, is shorter than the Ambartsumian- 
Spitzer time scale, f^^treiax, where fn > fN denote the fractional changes in total 
energy and number of stars per relaxation time, t^eiax- The resulting evaporative 
lifetime t^v — 'iO.btreiax for isolated clusters is consistent with Fokker-Planck and 
N-body simulations, where treiax is expressed in terms of the half-mass radius. 
We calculate the grey body factor by averaging over the anisotropic perturbation 
of the potential barrier across the tidal sphere, and derive the tidal sensitivity 
dlntev/dy ^ —1.9 to —0.7 as a function of the ratio y of the virial-to-tidal radius. 
Relaxation limited evaporation applies to the majority of globular clusters of the 
Milky Way with N = 10^ — 10^ that are in a pre-coUapse phase. It drives 
streams of stars into the tidal field with a mean kinetic energy of 0.71 relative to 
temperature of the cluster. Their S shape morphology leads in sub-orbital and 
a trails in super-orbital streams separated by SAa/Q in the radial direction of 
the orbit, where Q denotes the orbital angular velocity and a the stellar velocity 
dispersion in the cluster. These correlations may be tested by advanced wide 
field photometry and spectroscopy. 



1. Introduction 



Globular clusters re present some of the oldest stellar systems as indicated by their stellar 
ages of at least 12 Gyr ( jPeebles fc DickdllOGSt IVandenBerg et al.lll996l : iHansen et al.l 120021) 
consi s tent with the ages of white dwarfs in the thin disk of the Milky Wav (INoh &: Scalo 



19921 : ISalaris et al.ll2000l : iDeGennaro et al.l 120081 : llsern et al.l 120081 : iRenedo et al.ll2010l ). the 
oldest halo field stars such as HE 2523-0901 of 13.2 Gyr±2.7 Gyr (13.4±1.8 Gyr based on 
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U/r chronometry al one, iFrebel et al.l ( 120071 )) and the age of the Universe of 13.75 ±0.11 Gyr 
(jjarosik et al.ll2010l ). The population of clusters of the Milky Way as observed today is but 



a shadow of a primord ial population (jChernoff fc Weinberg Il990l : iMcLaughlin &: Fall 



2008 



Prieto fc Gnedin 2008 ). that may have formed in the era of reionization (jVandenBerg et al. 



1996[ ). Their ages hereby exceed their relaxation times by a factor of a few to an order of 
magn itude. They are essentially virialized close to thermal equilibrium (IMevlan fc Heggie 
1993) and some 80% satisfy the King model fKM. lKinej ffl966[):IChernoff fc WeinbereJ Jl99ohl. 
which bears out well also in numerical simulations ( IChatterjee et al.ll2010l ). The exception is 
a minor fraction of clusters in the direction of the galactic halo with a ceritrally peaked core. 



commonly designated as p ost core-collapse (PCC) flGhernoff fc Djorgovksilll989l : lTrager et al. 
I993I : Ijohnson et al.l[2008h . 



Prior to core collapse, the evolution of stellar cl usters is well-described by two-body inter- 



actio ns in the small angle scattering approximation (jSpurzem fc Aarsetlj|l996l:lAarseth fc Heggie 



1998 ). that naturally gives rise to evaporation of the cluster (lAmbartsumianI Il938l : ISpitzei 
19401 ). This proces s is amenable to detailed m odeling by Fokker-Planck and N-body simula- 
tions (jHenorull96ll : ISpurzem fc Aarsethll996r). there results of which can be interpreted sta- 
tistically (reviewed in lMeylan fc Heggid ( 1l997l )) in terms of evaporation, tidal effects, stellar 



mass loss, mass segregation and heating and crossings of the galactic disk ( Gnedin &: Ostrikei 



1997; 


Watters et al. 


2000; 


Josh et al. 


2001; 


Fregeau et al. 


2002; 


Heggie 


2001 


Chatter iee et al. 


2010; 


Heggie 


2011 


)• 



There is an increasing interest in these various processes to develop efficient semi- 
analytic approaches to the birth and e volution of large numb ers of stellar clusters as tracers 
of cosmological evolution of galaxies (IKruijssen et al.ll201ll ). Evaporation stands out as a 
common process in many-body systern s that includes for instance atmos p heric science (e.g. 
Charnberlairu (ll963r);lHavs fc Liul (11965 ) : IChamberlain fc Campbelll (119671 ) ; IShizgal fc Blackmore 
Jl986h : I Johnson et al.l J2008h : IjonnsoJ J201ohl. where i t extracts th e lightest elements from 



an otherwise nearly thermal state (e.g. IChavanis et al.l ( 119961 . |2002| )) 



Evaporation is a radiation process that defines a finite Kelvin- Helmholtz time scale due 
to energy extraction by the escaping stars. Key to their evolution, therefore, is an accurate 
determination of the energy loss carried off by stars escaping the cluster a determination of 
the time scale of evapo ration, in par ticular whether it follows the relax ation time scale o r 
the crossing time scale (jHeggidl200ll ) or possibly a combination thereof (IBaumgardtl I2OOII ) . 
The ener gy extraction by esca ping stars is potentially useful in classifying post core-collapse 
clusters (IKiipper et al.ll2008bl ). 



Evaporation of stellar clusters has received considerable attention (reviewed by I Jo. 



mstone 



(I1993I )) with an emphasis on mass loss typically more so than energy loss (but see 



Wielen 
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( ll967l )). Yet, the energy distribution function of all stars, those remaining in the globular 
cluster and those expelled by evaporation, continuously spreads towards increasingly neg- 
ative and positive values associated with contraction and heating. The kinetic energies of 
evaporated stars in the halo hereby carry a thermal imprint of the history of the clusters. 
Since the energy distribution function is unbounded in the point particle limit, there is the 
potential for extended lifetimes already at small binary mass fractions. However, observa- 
tional evidence for a dynamically significant population of compact binaries has not been 
demonstrated, nor in the correlation between th e observed populat ion of luminous X-ray bi- 
naries and the binary star-star interaction rate (jPooley et al.ll2003l ) that points to relatively 
small binary populations that do not necessarily affect the lifetime of the clusters. 

In this paper, we explore the parameter range for relaxation-limited evaporation, the 
effect of energy loss, and some observational consequences for the present state of globular 
clusters in the Milky Way. We discuss the thermodynamics of evaporating globular clusters 
with tidal perturbation in the large-N limit of equal mass stars as a function of a tidal 
parameter 



R 

y = —, 
n 



(1) 



representing the ratio between the virial radius R (close to the half-mass radius, to be 
discussed further in §2 below) and the tidal radius rj of the cluster. For old clusters such 
as of the Milky Way, rt is much larger than the radius of the cluster. The tidal sphere 
has the effect of lowering the potential barrier for stars to escape along orbits into the 
the halo of the galaxy and, generally, introduces stron g anisotropy in the directions of the 
escaping stars defined by the Lagrange points Li and L2 ( lHaylilll967l ). Existing approaches to 



incorporate tidal effe c ts commonly use the minirnum energy for escape a t Li (jSpitzer 



1987 



Lee fc Ostrikei] 119871 : IChernoff fc Weinberg! Il990l : iTakahashi et all 119971 : IRoss et all Il997h . 
effectively applying a spherically symmetric extension thereof to the tid al sphere. However , 
since last scatterings arise typically from within the half- mass radius (IBaumgardtl |2002| ) , 
the escape rate is effectively determined by the probability of escape derived from spherical 
averaging over the anisotropic tidal field for randomly oriented ejections. We thus derive the 
evaporative lifetime up to the moment of core-collapse to be 



te. ~ 20.53 (1 - 1.90y + 2.78y2 - 2A0y^ + 0.8%^ + 0(y^)) Ueia., 



(2) 



where the relaxation time tj-eiax expressed in terms of the half-mass radius Rh (the coefficient 
is 13.85 when treia x is ex pressed in R). While our approach is a largely explorative extension 
to lAmbartsumianI ( 119381 ) focused on the sensitivity of the evaporative lifetime to ener gy loss, 



some c onfidence can be gaine d by comparing (|2]) with known results review' ed in ISpitzer 



(119871 ): ITakahashi et al.l fll997D : iMevlan fc Heggid ( 119971 ): iHeggie fc HutI f l2003h . 
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Early in their evolution over the first few Gyr, globular c l usters may have experienced 
considerable mass loss in stell ar winds ([Aarseth fc Heggielll998l : iBaumgardt &: Makindl2003[ ) 



due to stellar evolution (e.g. lApplegatd (Il986[ )). that created a crucial facto r in the hkeli- 



hood of survival of a cluster after its formation (jchernoff fc Weinberg 1990 ). A common 
approximation used in Fokker-Planck and N-body simulations is an exponential evo lution 
on a time scale tm (lAarsethl Il985l : IChernoff fc Weinberg] Il99d : iMeylan fc Heggid 119971 ) . The 
same approach is readily incorporated in a thermodynamic model as a regular perturbation 
provided that » treiax- Evaporation has also been found to narrow the mass spectrum 
due to a relatively enhanced rate of escape of the li ghtest stars, by which older clusters tend 
evolves similarly to a system of equal mass stars ( IChernoff fc Weinberglll990l ). The afore- 
mentioned analogy with atmospheric science is useful here, noting that evaporation leads to 
an early escape of light molecules such as hydrogen and helium from the terrestrial atmo- 
sphere leaving only the relatively heavy molecules of nitrogen and oxygen to survive. Early 
evolution of star clusters therefore deviates from the evaporation driven evolution of single 
mass clusters. However, for the old clusters of the Milky Way at the present epoch, these 
considerations are less critical, and will not be considered here further. 

As an application, we focus on high velocity stars in an d around glo bular clusters that 
may provide "real-time" snapshot of last scattering events ( lHenonlll96ll ). The se last scat- 
tering events from within the half-mass radius (in the limit of isolated clusters, IBaumgardt 
( 120021 ) ) and most frequently deri ve from the co r e, where the den sity is highest and the local 
relaxation time is the shortest (lAarsethl Il974l : iJohnstond Il993l ). This includes tidal tails 
produced by ejectio ns that, indirectly , serves as a diagnostic for the c luster temperature as, 
e.g., in NGC 6254 Jheon et al.lboooh and 47 Tuc Jhane et allboiOal lbf). where the velocity 
dispersion rises to about 60% at the tidal radius relative to the velocity dispersion inside the 
cluster. Observational data on tidal tails hereby provides a signature th at may be used to caL 



ibrate numerical mo dels when ap plied to individual cluste r s, e.g.. M15 (|Grillniair et al. 



Drukier et al. 19981), Palomar 5 (Odenkirchen et al. 2001 : Mevlan et al.l 



by A^— body simulations (e.g. iHeggiel ( )2001 ): Lee et al. (j2006 )). This may include further 
processes such as additional heating by past disk crossings for clusters near the center of th e 



2001 



1995 



Dehnen et al. 



2004h . NGC 5466 dBelokurov et al.lbooel) and Palomar 14 JSoUima et al.lboill ) and their tails 



galaxy, but possibly less so for the outer clusters ( IGnedin fc Ostriker 



Some of these tida l tails display an 5— shaped morphology ( Heggiel 



1997 



2001 



Leon et al 



200 



Belokurov et al. 



2006 



Kroup a'ioos') along the Lagrange points Li and L2 on short scales, e.g., in NGC 5904 



( iLeon et al.li2000i: iJordi fc Grebell 120101) . and along the orbit of the cluster on long scales, 
e.g., NGC 5644 ( jPellhauer et al.ll2007l ). where the leading (trailing) tail is inside (outside) the 
cluster orbit around the galactic center. The temperature of the clusters v aries broadly, from 
about 1 km/s in Pal 12 (ILeon et al.ll2000l ) to about 15 km/s in a;-Centauri (Ivan der Yen et al. 



- 5 - 



20061 ). with evidently different prospects for tlie manifestation of tliese tidal tails. 



The planned Large Synoptic Survey Telescope (ILSSTI 120091). the Baryon Osci 



lation 



Schlegel et al.l fl2009f)1 th e 



Spectroscopic Survey (BigBOSS, a multi-object spectrograph, 
Japan Astr ometry Satellite Missio n for INfrared Exploration (JASMINE lGouda et al.l ( l2010l )). 
and GAIA (ILindegren et al.lll995l ) seem to be ideally suited for studying these tidal tails by 
large field of view photometry and spectroscopy. This development offers novel opportunities 
for testing the physics of evaporation by measuring correlations between the kinematics of 
stars in tails to the thermal state of the clusters. 

In §2, we define notation and the basic quantities describing the thermodynamic state 
of virialized star clusters, their crossing and relaxation time. In §3, we introduce the gen- 
eral set-up for the evolution equations, in §4, the parameter range for relaxation-limited 
evaporation is identified. In §5, we calculate the evaporation coefficients for the number 
of stars and energy emitted per unit of relaxation time. In §6, we extend our model by 
regular perturbation analysis to include tidal interactions. In §7, we discuss the remaining 
relaxation-limited evaporative lifetime of the globular clusters in the Milky Way. We give 
model predictions for observational test on tidal tails in §8 and summarize our results in §9. 



2. Virialized star clusters 



By the large number N = 10^ — 6 of stars, the globular cluste rs of the Milky Way ar e 
essentially virialized with constant mass-to-light ratio in their cores ( iMeylan fc Heggielll997l ). 
For globular cluster s of the Milky Way , corrections to the Virial Theorem due to tidal forces 
are generally small (ILee fc Hoodlll969l ). In light of the aforementioned KM morphology, we 
begin with a homogeneous and isolated cluster of single stars of the same mass, m, and no 
binaries. The system is then described by a total number of stars, and virial radius, R. 
The Newtonian binding energies Uij = —Gm'^/\ri — rj\ between stars i and j (i ^ j) give 
rise to a total energy H representing the sum of total potential energy U between the stars 
and the total kinetic energy of their individual motions, satisfying 



H = U + Ek 



4 



where we used the Virial Theorem, 2Ek + U 



^ 

= 0, on the mean values 



U 



2N 



GmM 
2R '■ 



Eh 



N~^Ek 



(3) 



(4) 



and where ([3]) is a defining relationship for the virial radius R (cf. iMeylan fc Heggid ( 119971 )). 
The ratio M/ R is a hereby a measure for the mean kinetic energy in the system. 
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The relaxation time for a homogeneous population of stars arises out of elastic 2-body 
encounters, produce random deflection angles Acf) ^ ±Uij/Ek, where Uij = —Grn?/h denotes 
the binding energy associated with the impact parameter h for an interaction between stars i 
and j. High velocity encounters hereby experience relatively small angular deflections, giving 
rise to asymptotic freedom generic to scattering problems. In a crossing time 

^3/2 



tc. 



R 

V 



up to a factor of unity, the numerous encounters of a star wandering through a cluster give 
a non-zero expectation value of the variance of these deflection angles. (The crossing time 
corresponds to the free-fall time scale within a factor of t wo.) The resu lt can be summarized 
in a relaxation time for an expectation of unit variance (jSpitzerlll940l ) 



t 



N 



relax 



-X 



2.4 X . ,3 Myr, 



(6) 



61nA^ " hiN 

where R = Ri pc. The approximations involved allow for slight variations, while preserving 
the main dependencies on the overall characteristics of the stellar system. The details of 
these differences shall not concern us here. 

The virial radius R is closely relate d to the half-mass radius are closely relate d, e.g., 
based on the Plummer model, we have (lAarseth fc Falllll980l : iMeylan fc Heggid 119971 ) 



2R 



-U 



0.38- 



(7) 



whereby R = 1.30rh. Consequently, trh 
estimated on the basis of rt- 



32N{rh/1.30) 

0.67treiax, where trh denotes the relaxation time 



3. Equations for evaporative loss of mass and energy 

Evaporation is a radiation process in which stars with sufficient energy escape, out to 
infinity for isolated clusters or out into the tidal field of the host galaxy. This process can be 
described by fractional changes and fn in the number of stars, N, and, respectively, the 
energy of the cluster per unit of a characteristic time, that is generally close to the relaxation 
time. It represents the statistical outcome of a relatively large phase space of possible stellar 
trajectories, that emerges from, most frequently, soft two-body interactions. 

When the relaxation time is intermediate between a dynamical crossing time and the 
evaporative lifetime of the cluster, the escape probabilities can be calculated in a thermody- 
namic limit subject to the Virial Theorem and radiation boundary conditions. In the absence 
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of thermodynamic equilibrium, the ratio Jh/Jn is anywhere from zero, e.g., when all stars 
are trapped or escape with at most the escape velocity, to arbitrarily large, e.g., in the un- 
likely event of a star escaping with a kinetic energy on the order of the total gravitational 
binding energy of the cluster. 



AmbartsumianI (119381 ) described evaporation of isolated clusters by a singular perturba- 



tion away from thermal equilibrium given by the tail of a Boltzmann velocity distribution 
representing a flux of stars with positive energy to escape to infinity on the relaxation time 
scale ^re/ax? that is. 



dN 
It 



where /at — 1/135. (Here, /a? is sometimes referred to as ^e, e.g.. ISpitzerl (Il987| )). Thus, 
defines the Ambartsumian-Spitzer time scale 

tAS ~ In t relax- 



(9) 



The true evaporative lifetime of cluster is shorter than ([9]), however, by continuous heating 
of the cluster in a fully nonlinear evolution, upon including the energy radiated away by the 
escaping stars. In addition to treiax, the characteristic time scales for a bound cluster, there- 
fore, may further include the crossing times tcr,i = {R/vt, Rh/vt,rt/vt}, where Vt represents 
the mean thermal velocity of the stars in the cluster. These considerations generalize (IH]) to 



— — fnN F {treiax jtc 



with accompanying loss of total energy 



dH 



relax J ^crA 



(10) 



Early attempts to include energy loss use order of magnitude scaling arguments (ISpitzer 



19401) or bv explici t calculation of energy frac tions in the tail of the velocity distribution 
(jvon Hoerner 1958 : King 1958 : Ossipkov 2008). This approach seems justified in the KM 
morphology, but requires careful estimation of the numerically small fractions involved, that 
repre s ent a singu lar perturbat i on a^v ay from thermal equilibrium. The early estimates of 
King! (119581 ) and Ivon Hoernerl (119581 ) give lifetimes shorter than 135treiax, but still much 
larger than what is found in numerical simulations by a few hundred percent due to an 
underestimate fn < /n- Present numerical estimates of the time to core-collapse range 
from tev = 10 — 20trh in terms of the relaxation time associated with the half-mass ra- 



dius, Tfi, that are considerably smaller than the lAmbartsumianI (119381 ) evaporation time, by 
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Fokker Planck and N-body simulations ( 


SDitzer 


& Shapiro 


1972; 


Cohn 


1980; 


Inaeaki 


1985 


Lee & Ostriker 


1987 




Spurzem & Aarseth 


1996; 


Drukier et al. 


I999I: I 


ieeeie & Hut 


2003) 



Here, rel atively longer t imes are found in King models with medium a nd low central concen- 
tratio ns ( lQuinlanlll996[ ) and intermediate times in a Plummer model (jChernoff fc Weinberg 
990), while short times are obtained in multi-mass models with broad mass distributions 



Chernoff fc WeinbereJ E990 : |Johnstonelll993l ) with minor dependency on velocity anisotropics 



(ITakahashil 119951 ) ■ 



In this paper, we shall formulate the evaporative lifetime in terms of the coefficients 
(/at, fn), derive ([2]) and determine its sensitivity to energy loss, going beyond that predicted 
by alone. To pursue this, we first identify a critical number of stars, for evaporation to 
be limited by relaxation. 



4. Relaxation-limited evaporation when the number of stars is large 

At a given location in the cluster, the Boltzmann distribution n{E) = N'(3e~^^ for the 
kinetic energy E (equal to the total energy as measured by a local observer) is the limit 
for detailed balance in a local thermal state, where A^' = n{E)dE denotes the total 
number of particles in a local neighborhood, f3 = l/ksT the inverse temperature and ks 
Boltzmann's constant. Relaxation of small perturbations away from the Boltzmann takes 
place on the relaxation time scale treiax associated with small angle scattering in our present 
approximation. The normalized distribution n = e^^n evolves to a uniform distribution in 
E. The process of relaxation for n is therefore equivalent to a diffusion process in a medium 
with uniform affinity (a vanishing chemical potential throughout), whereby dtu = t~^iax^E^y 
i.e., 

dME,t) = t;^l^{dln-PdEn). (12) 



Evaporation is a process in which particles of energy E > E^, where Ec depends on the 
location in the cluster, escape on the time scale of the crossing time, tcr, set by the size of 
the system, that provides an additional relaxation time for radiation leaking out to infinity. 
With the scaling x = E/Ec, a = Ec/ksT, 

ftn(E)dE' 
we are led to consider the diffusion equation 

9tn = t-^iax [din - ad^n) - t'^Uix - l)n, (14) 
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where U denotes tthe Heaviside function. About a thermal distribution, Generally, the 
evolution of n = n{x,t) depends on both timescales itreiax^icr)- The case of treiax » tcr 
represents slow relaxation, relative to which evaporation on the time scale tcr is essentially 
immediate. It will be appreciated that f lT^ is similar in form to an inhomogeneous Kompa- 
neets equation. The bound for this asymptotic result to hold allows by balancing the second 
and third terms in ([HI), i.e., 



trel 



ax 2 

~ a 



te 

with a corresponding critical value N = Nc satisfying Nc — 2181ogA''c, i.e. 



(15) 



1600 



(16) 



by ([5]), (|6]) and the lAmbartsumianI ( 119381 ) value a ^ 6. Thus, evaporation of clusters with 
N » Nc is limited by relaxation, which justifies the starting point dHHn]) and hence (llOtilip 
with 



t 



relax 



(17) 



and the corresponding Kelvin-Helmholtz time scale 



relax • 



(1^ 



The parameter range for relaxation-limited evaporation ( II 7p is amply satisfied by the globular 
clusters of the Milky Way with = 10"^"^ Young open star clusters in the disk of the Milky 
Way can be different with A^ as low as a few hundred. For these systems, a full equation of 
the type (fT4|) must be used. 

A gradual departure away from the linear scaling of ev aporative lifetimes oc treiax has 
been observed in N-body simulations by iBaumgardtl (120011 ) in the range A^ = 10^ — 10^, 
which we here identify with the lo w-N regime o f (fH|) defined by ( 1T6|) . where evaporative 



lifetimes depend on treiax and t cross- IBaumgardtl ( 120011 ) models the low-N scaling with 



relax cr 



6\nN 

N 



t 



relax i 



(19) 



where /3 ~ 0.25. However, Eqn.(12) in IBaumgardtl (l200ll ) from which (fT9l) derives, does 
not recover the Boltzmann distribution in the absence of evaporation (defined by tcr oo). 
Hence, it fails to predict a corresponding critical value in A^ that explaiii s this t ransition. By 



([T6l). w e a re able to ascertain that the discussions of IBaumgardtl (120011 ) (and iLamers et al. 



(l2005al lbl): iGieles fc Baumgardtl ( l2008l )) pertains strictly to low-A^ asymptotic behavior in 
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N < Nc with no bearing on high-A^ asymptotic behavior relevant to the present-day globular 
clusters with KM morphology in the Milky Wa y (with most having relatively small tidal 
parameters ([T])). A similar conclusion is given by iMcLaughlin fc Fall! ( 120081 ) . 



Baumgardt 


( 


2001 


); 



The above determines the transition to be around (fT6!) for A^. The precise value of is not 
cr itical, as the transition is rather smooth, as may be seen from numerical results in Fig. 4 



of lBaumgardtl (120011 ). 



In what follows, we shall continue in the parameter range N » Nc, where the evapo- 
ration process is governed essentially by tj-eiax- 



5. Evaporation coefficients and fn of a Maxwell distribution 



For an isolated cluster, an individual star is bound to escape whenever its energy is 
positive, i.e., ek + u> 0, where denotes the kinetic energy and u the gravitational binding 
energy to all other stars in the cluster. For star i, we have u = J^jLiUij (sum not including 
j = i). B y (HI) and the Virial T heorem, the mean kinetic energy of the escaping stars hereby 
satisfies (lAmbartsumianl Il938l ) 



-u 



-2U = AEk. 



(20) 



Correspondirigly, th e escape velocity is about twice the mean velocity of the stars (lAmbartsumian 
19381 : ISpitzeiill940l ). For a thermal distribution, it represents a fraction 



N 



Jo 



,3 o2 



e 2' 



{A7Ts^)ds 134.44 



(21) 



described by a Boltzmann distribution, where = Ek/E^, E}^ = ^ksT, where T denotes the 
temperature. Alternatively, we may identify the cluster with the collection of bound stars 
whose velocities are restricted to the truncated Boltzmann distribution. In this event. 



N 



/■oo 
2 e 



1 



(22) 



where k = 1.452165 represents the modified relation E/^ = nkBT obtained from ( 1T3|) . i.e. 



r 

Jo 



2 „2 



s e 



(47rs2)rfs 



/o e~'^'*^(47rs2)(is 



1. 



(23) 
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The difference in fl2m22p may be considered representative for the uncertainty in the present 
thermodynamic approach, and we will proceed with both to keep track of uncertainties in 
any of the derived quantities. 

The relaxation time ([6]) is the instantaneous time scale for the relaxation of a stellar 
system, given its current state. It defines the rate at which the high velocity tail in the 
Boltzmann distribution of escaping stars is replenished by upscattering of low velocity stars, 
giving rise to the evaporation rate ( llUp and the instantaneous evaporation time or 
21% less based on (1221) . The estimate represents an extrapolation based on the initial 
conditions of the system, which otherwise evolves nonlinearly over multiple relaxation times. 
Because stellar systems are self-gravitating, their temperatures increase as they loose stars 
by evaporation. Therefore, gives an upper bound for complete evaporation, not the true 
time to complete evaporation in a fully nonlinear evolution. 

In the process of evaporation, the system evolves also in response to energy loss due to 
positive energy of the escaping stars. A complete system of equations, therefore, includes 
energy balance. Analogous to previous arguments, the total rate of energy loss ('luminosity') 
changes the total energy H of the cluster at a rate (ITT]) with 

_ /r.^e-F(W)d. _ 1 

s^e--2^\A7rs^)ds 28.75' ^ ' 

Here, we identify H with the total energy of the stars below the threshold (all the stars that 
make up the cluster) and fn with the fraction of total energy carried off by the evaporating 
stars, that replenish the tail above the threshold in each relaxation time period. 

Note that fn is calculated by the relative fraction in kinetic energy in the tail of the 
velocity distribution, representing a flow in momentum space from local interactions (in 
coordinate space) between "low" velocity stars (below the threshold, in the cluster), where 
kinetic energy and potential energy are virialized. Thus, replenishing the tail above the 
threshold with kinetic energy is tightly correlated to H defined by the stars below the 
threshold and vice versa. Following (!22|) . a more precise definition is, therefore, 

/o%2e— '(47rs2)tis 23.74' 

which differs from (!24|) by 17%. 

The system (llOtllip is closed by ([3]). Integration of the two ordinary differential equations 
gives 



(26) 
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^^ fN + fn ^ 5.71, a = ^ + 1 = 6.71 (27) 

JN 



where 



based on 7 = 5.712 from f l2T|) and (1241) . As a ratio, the estimate fl27j) is stable in the first 
two digits in view of 7 = 5.736 as follows from ( l22l) and ( l25l) . The correlations ( 12611271) are 
independent of the choice of evolution time scale F. 

As illustrated by ( 12T]) . fl22|) . ( |24|) and f l25|l . evaporation is a singular perturbation away 
from thermal equilibrium which introduces high sensitivity (of the coefficients {fNifu)) to 
the critical escape velocity, even though some dimensionless results are relatively stable such 
as (I27p . For example, a variation of +1% in the threshold defining the tail gives rise to a 
variation of about +10% in tev f lM|) . For small changes A, we have 

, /2°!Ae-i^'(47rs2)rfs 1 + 5.57A + 13.92A2 + 0(A3) 



_ jZ^s\-'^^\A7:s^)ds _ 1 + 4.73A + 9.45A^ + 0(A3) 
^""^ ^~ JZ s^e-l^\A7rs^)ds " 28.75 ' ^''^ 

Simular results hold for the fractions (f^Jn), i-e., /i^//jv(0) ~ 1 + 5.3781A + 12.9308A2, 
fn/fHW ^ 1 + 4.5419A + 8.6492A2. Fig. □ shows both (/^, fn) and (/;^, /^) as a function 
of A. Corresponding to ( I28ll29l) . the evaporation time satisfies 



re,(A) = 12.82 (1 + ) (1 - 5.01A + 14.13A2 + 0{A^)) (30) 

V In A^o / 

and, similarly, re„(A) = 10.60 (l + (1-4.82A + 13.15A2+0(A3)) in the approximation 

defined by {f^,fH)- Consequently, we have 

^ = 4.71 (1 - 0.8434A + 0.23094A2 + 0{A^)) (31) 

JN 

and similarly §- = 4.74 (1 - 0.8460A + 0.2150A2 + ©(A^)), which evidently is essentially 
the same for (/tv, fn) and (/^, /^). We conclude that the impact of evaporation on energy 
is about five times larger than that on mass. 



5.1. The evaporative lifetime of an isolated cluster 



For an isolated cluster, integration of (llOtilip with ([H]), fl2B]) and F = t^^^^^ gives 

d^ / In/i 
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Fig. 1. — {First and second column.) Shown are the evaporation fractions {fNifu) and 
Un) fh) a function of the relative departure A from the nominal escape velocity {circles) 
over the range A = [— i, i], including quadratic approximations thereto {continuous lines), 
which differ by about 20%. {Third column.) The resulting 7 = 1 + -^ ( 1271) are essentially 
the same for both pairs of fractions {continuous and thick lines, respectively). The normal- 
ized evaporation time in shows a steep dependency on A due to the faster than 
exponential decay in the tail of the Boltzmann distribution. 



1 3 

with r = t/t°eZaa;! whcrc t°eZax = "^kiM)" '^^'^'^^^^ initial relaxation time defined by A^o and 
the initial radius Rq, /i = the normalized particle number and g = 2 + |7 = 10.57. A 
finite-time singularity solution obtains by ignoring the logarithmic term in ( l32l) . 



N = No[l 



T 

T ev 



(33) 



where t^v — Q'^f^^ = 12.82. A next order approximation for obtains by including a 
time average of the logarithmic term following (133|) . With In/i ~ —tJn for most of the time 
(away from the singularity), the time to complete evaporation extends to 



N 



1 



2q In A^o 



12.82 1 



0.05 
InA^n 



(34) 



For all practical purposes, the evaporation time is about 13 times the instantaneous relax- 
ation time, apart from a weak positive dependence on the initial particle number A^o by no 
more than 5%. 

Solutions of the type (l33l) are canonical for a broad range of radiation processes with 
generally diff erent coefficients, that would include evaporation of relatively low N star clusters 
according to iLamers et al. J2005bh on the basis of (1191) . 
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The evaporative lifetime (l34j) serves as a lower bound to the core-collapse time of iso- 
lated clusters. Quite typically, scattering of stars occasionally has an impact parameter much 
smaller than bmm = Gm/v"^, in which case the interaction can be inelastic by tidal interac- 
tions that may lead to tidal capture and binaries formation (e.g. ( iPooley et al.ll2003[ )). This 
dissipative effect is ignored in idealized point-particle N-body simulations. While numerical 
simulations suggest that binary fractions of a 0.1% can alr eady produce an order of magni- 
tude increase in the lifetime of the cluster (I Josh et al.ll200ll ). a detailed study identifying the 
threshold for a population of compact binaries to affect the lifetime of a cluster appears to 
be lacking. 



6. Stimulated evaporation by tidal fields 



In clusters bound to 
gravitational field of the 
Lagrange points Li and 
depends on the direction 
deflection due to the angular m omentum of the g. 
N-body simulations (Fig. 
given by the tidal radius (ISpitzer 



a host galaxy, the escape velocity will be perturbed by the tidal 
galaxy, which generally facilitates stars to spill over across the 
L2, i.e., Cfc + M > —6 following (120]) . where 6 > 0. Ah initio, 6 
relative to the line passing through Li and L2 with an additional 



1987 



jular cluster, as can be seen in direct 



7.4 in iMevlan fc Heggid (|l997l)'). The distance to Li aiid L9 is 



Chernoff fc Weinberg! Il990l : iTakahashi et al.l 119971 ) 



rt 



Mr 



3M 



G 



(35) 



where Mq denotes the mass of the host galaxy and a the separation of the cluster-galaxy 
system. 

Existing approaches to incorporate tidal effects commonly use the minimum energy for 



escap eatl^iJlSmtzer 



9871 : iLee fc Ostrikei]ll987l : IChernoff fc Weinberdll990l : ITakahashi et al. 



19971 : IRoss et al.lll997l ). effectively applying a spherically symmetric extension thereof to the 
tidal sphere. However, the probability for escape depends on the orientation relative to the 
anisotropic tidal field, and hence requires averaging over the directions of last scattering 
events, here illustrated by F igs. [3lll] below. The possibility of long hved delays in escape 



(IBaumgardt fc Makindl2003l ) conceivably modifies the escape probability, but this effect will 



be small in the limit >> A^^^. 

As the relaxation time is smallest in the high density central region of a cl uster, most 
of th e escaping stars derive from high velocity jumps produced in the center ( iJohnstone 
I993I ). The tidal perturbation of the evaporation process can hereby be effectively modeled 
by studying the escape of stars that are emitted from a region about the center that is small 



- 15 - 



relative to the tidal radius of the cluster (relative to the distance be tween the L a grang e 
points Li and L^) in essenti ally randomly oriented directions; see also iBaumgardtl ( 120021 ): 
Baumgardt fc Makinol J2003h . 



The trajectory of a single star in the potential outside of the core of a globular cluster of 
mass Mc = NM is the solution of a restricted three-body solution described by an effective 
potential in a Cartesian frame of reference {x,y,z) = (oc + r cos^,r sin^cos0, rsin^sin0) 
with angular velocity Q,^ = '^^g+'^c ^bout the 2;— axis, given by 



d, 



1 /etc 
1 + - — 

2 V a 



(36) 



where a = y^x^ + is the distance to the axis of rotation, and = ■\/(a; — a^) 
and Tg = iy(x~+~a^p~+~^M-^ denote the distances of the star to the center of mass of 



the cluster and, respectively, the host galaxy. The latter have distances 



a 



Mg + Mc 



and 



-a- 



to the center of mass of the cluster-galaxy system, here at the origin of 

illustrate the tidal 



(x, z). We use units in which Newton's constant is equal to 1. Figs 
perturbations about the critical escape velocity, showing the critical onset of evaporation 
from the tidal sphere as a function of stellar velocity and anisotropy produced by the tidal 
field. There is further a gradual trend towards reduced anisotropy with increasing velocity. 

The tidal interaction generally lowers the barrier for stars to escape from the tidal sphere 
by an energy perturbation 5 that depends on the direction, as 5(r, 9, (p) has period vr in 6* 
(27r in 0) with a pronounced amplitude in 6 and minor variations in shown in Fig. O The 
Virial Theorem for the cluster fl20l) now becomes 



2Ek + <l> = 0, 



(37) 



where ^ = U — 6 with U and 6 defined in (jl]), respectively, fl5^ . It follows that U = —2Ek + S. 
The condition on the kinetic energy Cfc of a star to escape becomes Ck + u = ek + 2U — 6>0, 
giving a critical kinetic energy AEk ^1 — 4^ j reduced by the ratio of S relative to the kinetic 

energy AEk required for escape if the cluster were isolated. As the core of a globular cluster 
is much smaller than r^, the escape of high velocity stars from it is effectively isotropic and 
the net escape rate obtains from averaging over all directions. The average of 6 over the 
tidal sphere satisfies 



1 

An 



6{rt,e,(f)) smeded(j) = r 



'0 JO '^t 
where the grey-body factor of the evaporation process satisfies 



(38) 



r ~ 0.11 - 0.12 



(39) 
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Fig. 2. — {Top.) Formation of S*— shaped streams shown at two times Qt = 10 and Qt = 40 
in the corotating frame of the cluster by the trajectories (blue) of 20 stars ejected uniformly 
in all directions from the cluster with mass M = 10~^Mg' with orbital angular velocity Q 
around the host galaxy of mass Mq- The trajectories emanate from the center at tq = O.lr^, 
where is the tidal radius = 1.5 {red; 2rt in green), relative to a distance of 100 from 
the center of mass of the cluster-galaxy system located at the origin. The results are shown 
for an initial velocity vo which is 0.943 times the escape velocity vq of the same cluster 
when viewed in isolation with no tidal field. {Bottom.) Transformed to the rest frame of the 
galactic center, the streams form extended tails as they would aggregate into diffuse rings 
in the absence of destruction by crossings with the galactic disk on the orbital time scale of 
the cluster {not shown). 
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Fig. 3. — Shown are the trajectories [blue) of 30 stars ejected uniformly in all directions 
of a globular cluster similar to Fig. |2]with initial velocities vo increasing by 0.5% in each 
window, normalized to the critical velocity -i/— of the same cluster in isolation with no 
tidal field. The tidal sphere forms an effective threshold for stars to escape. The results 
serve to illustrate the critical onset of anisotropic evaporation of high velocity stars in the 
tail of an extended velocity distribution in the presence of a tidal field. 
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Fig. 4. — Shown are the trajectories as Fig. |3l now with a choice of initial velocities vq 
increasing by 2.2% in each window. The results serve to illustrate escape across the tidal 
sphere with an increasingly isotropic outflow of stars with increasing energies. 
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in the limit of small ratios M^Mg as shown in Fig. [51 With A = 2 — 2y 1 — and (jl]), 
we arrive at 

A ~ ^ ^ r-. (40) 

We are now in a position to calculate the modified evaporation rates due to a finite tidal 
radius r^. 




Fig. 5. — {Left.) Shown is the fraction of escaping stars from a cluster with mass M^/Mg = 
10~^ {6 X rt/M = 11.86%) and tidal radius r^/a = 1.5, where a denotes the to the host 
galaxy, computed according to Figs. [SHU The initial velocities v are normalized to the 
critical velocity for escape Vc if the cluster were isolated, whereby 1 — {v/vcY expresses the 
reduction in the potential energy barrier for escape due to the tidal field. The results show 
a reduced threshold in the kinetic energy for escape by about 14% and a reduction of less 
than 12% for a 50% probability of escape. This sensitivity of a 50% change in the fiux of 
stars escaping in response to a 2% change in A is consistent with the theoretical results 
summarized in Fig. [H (Right.) Shown is the orientation averaged < 6 >= 6 kinetic energy 
for escape from a cluster bound to a host galaxy, relative to the kinetic energy for escape 
if the cluster were isolated. For typical cluster-to-galaxy mass ratios M^/Mg « 1, the 
resulting grey-body factor F = 5 x rt/Mc asymptotes to about 11%. 



6.1. The evaporative lifetime of a cluster in the tidal field of its host galaxy 



For the globular clusters of the Milky Way at the present epoch, the tidal ratio [T] 
assume rather moderate values broadly distributed around 0.1 (IBaumgardt et al.l l2010l ). 
For this reason, we purse a regular perturbation of the evaporative lifetime in this tidal 
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ratio, as opposed to a fully nonlinear treatment that would include tidal ratios of order one 
dCieles fc Baumgardtll2008[ ). 



Following the previous section, we integrate the ordinary differential equations for mass 
and energy fllOtllip with closure (/at,///) given by fl28tl29p and A in ([1]) expressed by (HOj) . 
For (13T|) . we have, by the leading order solution (126|) . 



^ = 4.71 (1 - 0.097?/ + OmSly^ + 0{y^)) = a f 1 - byoi^^ + 0{y^ 
In ^ 



(41) 



where we take into account only the linear dependence on y, a = 4.71, b = 0.097, and yo 
refers to the initial value. Then the leading order solution for the energy evolution is 



H = Hon "e5+3« 



-4.71 



in 

' 3 



(42) 



where c = 0.072?/o, and similar to (!26|) we find 



N 
iVo 



-7 



1-C 



iVo 



R_ 

Ro 



N 



1 + c 



(43) 



the correlations fj43l) are independent of the choice of evolution time scale 



T 

n~~ 

Similar to ([SMS 
F. 

As when deriving fl52]) in the previous section with fi = N/Nq, the evaporation time 
relative to the initial relaxation time follows by integration from ( 1TU|) with (P) and, now, 
(143]) instead of (|26D, 



■ In/i 3c 

1 H ^- u 3 

IniVo 2^^ 



;i + OMyofi-^ + 1.62/o> 



(44) 



where 



1 

135 • 



Fig. Elshown the result of numerical integration of (jH]) using a canonical value Nq = 10^, 
where we infer a weak logarithmic dependence on Nq by comparing results with A^o = 50. 
With a polynomial fit to the tidal dependence, we arrive at ([2]), where we dropped the 
subscript to initial values of y and A^, and where we dropped a factor (1 + 0.002 In Ag)"^ 
with extremely weak dependence on A^ = A'g x 10^. We emphasize that ([6]) refers to finite- 
time singularity solutions of single-mass clusters, as a regular t idal p erturbation of isolated 



clusters in their pre-coUapse state (e.g. Fig. 1 in iBaumgardtl (120021 )). The detailed shape 
of the graph of /x(t) = N{t)/No in case multimass star clusters can be quite different with 
a continuously evolving mass function due to relatively fast evaporation of low-mass stars 
(e.g. Fig. 1 in Baumgardt fc Makino ( 20031 )). 
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The numerical result shows a tidal sensitivity 



d In Tp. 



dyo 



-1.9, 



d In r„, 



-0.70, 



(45) 



giving rise to reductions up to a factor of a few when ?/o approaches 1 as shown in Fig. El 
Very similar results follow from (/^, f^). For young clusters with possibly |/o > 4 such that 
A > 0.2, (144]) acquires contributions from A to higher order than those given in fl28tl29p . 
which would require a further refinement of numerical integration. 




y =0..1(N =10°) 
yj,=0..1(Ng=50) 




numerical integration (Njj=10 
polynomial fit 

numerical integration (Njj=50) 



4 I 1 \ \ 1 \ \ 1 \ \ 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fig. 6. — {Top.) Evolution of the number of stars /x = N/Nq for various initial tidal param- 
eters I/O according to (jS]), where r denotes the time in units of treiax defined by the virial 
radius R. (Bottom) The slope to r^^ is —0.94 at yo = 0, showing strong tidal sensitivity, 
here by numerical integration of ( l44l) . The results are essentially the same for different A^^o, 
showing that the process of relaxation-limited evaporation is essentially independent of N. 



6.2. Comparison with numerical simulations 



We may compare the above with existing approximations for th e evaporation rate of 
stars in tidally truncated clusters (ISpitzerl 119871 : iTakahashi et al.l 119971 ). that are focused on 
/at with point wise extensions of the potential barrier at Li, as opposed to spherical averaging 
( l38l) over the tidal sphere. This generally over-estimates the sensitivity to tidal interactions 
by a factor of about 4 given the amplitude of about 0.5 in ^6 as shown in [51 Furthermore, 
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5]) may be compared with that resulting from the fit 



6.67 



(46) 



to numerical s imulations in terms of the half-mass radius (jSpitzer fc Chevalieii Il973 
Auguilarlll988l ). satisfying 



din tp 



-1. 



(47) 



2/0=1 



Our result given by fl45|) shows a generally y— dependent tidal sensitivity of the evapo- 
ration lifetime, that is overall consistent with the fit fH7|) . The resulting evaporation lifetime 
(j2]) connects yo = I smoothly to the finite lifetime for isolated clusters {yo = 0). It will be 
noted that the fit (H6l) is restricted to a neighborhood of ?/o = 1, as it predicts an infinite 
lifetime for isolated clusters or, equivalently, clusters with very small tidal parameter ([1]). 



Gieles fc Baumgardtl (120081 ) discuss a polynomial approximation to the instantaneous 
evaporation rate as a function of tidal radius. Beyond their cut-off of y = 0.05, their Eqn.(7) 
predicts a logarithmic derivative of -1.5 for the evaporative lifetime to variations in the tidal 
ratio, similar to ( l45ll derived above. 



7. The remaining evaporative lifetimes of globular clusters of the Milky Way 



For globular clusters of the Milky Way, we recall that most have KM morphology, 
all satisfy (fT6ll and y << 1, by which ([2]) is applies with rather minor dependency on tida l 
effects. Independ e nt arg uments pointing to the same are given by lMcLaughlin fc Falll ( l2008l ): 
Prieto fc GnedinI to08( ]. 



Clearly, the tev of the globular clusters exceed a Hubble time in agreement with their 
ages, where tev = 14 Gyr separates the half-mass parameters on the right and the core- 
parameters {Nc = {Rc/ R)^N, Rc) on the left. Overall, the distribution follows a trend along 
the paths of steepest descent of treiax, where the core-collapsed clusters continue this trend 
well below the model limit N = 2. 

Fig. [7] also shows that there is no correlation between the X-ray luminous binaries 
mentioned in the introduction and globular clusters with or without prior core-collapse, 
except that clusters with these X-ray sources tend to be of high mass for a given cluste r radiu s. 



If (compact) primordial binaries contribute to a PCC morphology (e.g.. iGao et al.l ( 119911 )) 



which is unclear at present, they are unlikely to contribute to the X-ray sources. Instead, 
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Fig. 7. — Shown i s the mass-radius distribution of 119 clusters in the Milky Way (from 



Gnedin fc Ostrikerl ( 1l997l )) in terms of (A^, R), where N = M^/Mq and R denotes the radius 
at half mass Mh, differentiated in r aorphology for King Model ( KM, circles) versus Post 
Core-Collapse (PCC, squares) ffrom lchernoff fc Djorgovfeil Jl989k adding FCC to NGC106 
and NGC4499). We add the distribution of inferred core parameters (iVc, R^ in the Ansatz 
Nc = {Rc/R)^N (KM, >; FCC, *) which, for FCC, defines a lower bound on the associated 
Curves of constant t^v are shown for 35 Myr, 3.5 Gyr, 7 Gyr and 14 Gyr. The 14 Gyr 
line separates the (A^, R) from {Nc, Rc) without regard to morphology, consistent with the 
present age of the Universe and the age of the clusters. The average tev for the clusters 
cores is about 1% of that of the clusters. The combined distributions of (A^, R) and {Nc, Rc) 
follows the path of steepest descent {thick blue line, slope 2.64) in treiax (two curved black 
lines), where the core-collaps ed clusters continue thi s trend well below the N = 2 model 
limit. A list of X-ray sources (IVerbunt fc Lewinl l2005l ) is included {red diamonds), showing 
some preference towards high N for a given R. Four clusters are between the curves of 14 
and 28 Gyr, that will evaporate in the next Hubble time consistent with the evaporative 
destruction rate of lAguilar et al.l (119881 ). 
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the X-ray sources are likely produced in the course of stellar relaxation by tidal captures 
( jPooley et al.ll2003l : iBenacquistal l2006l ) with no need for a primordial binary population. 



8. Cluster-tail correlations produced by evaporation 

Evaporation of stars modulated by tidal fields give rise to anisotropic outflows with 
bipolar morphology, commonly referred to as tidal tails. Relaxation limited evaporation 
provides a launching mechanism that predicts a correlation between the velocity of growth 
of these tidal tails an d the v elocity dispersion in the cluster as defined by its temperature 
( iKiipper et al.ll2008al Jbl. 120091 ). We here quantify this correlation suitable for detailed obser- 
vational studies. 

Following ( ITOlfTT]) . the mean kinetic energy efc(r) of the escaping stars at distance r are 
related to the mean kinetic energy of the stars in the cluster by 



efc(r) 



— -A\Ek + —^y-: 

Jn J R r 



where we encounter the evaporation efficiency 



^ - 4 = 0.71 (l - 0.097?/ + 0.0031|/2 + 0{y^)) 
In 



(48) 



(49) 



by dUD. 



Clusters with y <! will produce a bipolar outflow of stars that is amenable to detailed 
observations as mentioned in the introduction. These tidal tails carry a record of the most 
recent evaporation process by the number and kinetic energies of their stars as measured 
relative to the center of mass of the cluster. The planned and upcoming LSST, BigBOSS, 
JASMINE and GAIA offer unique prospects for extracting kinematic information in a wide 
field of view around globular clusters with angular resolutions around 10 ^ as, resolving 
velocities of individual stars to about 1 km/s at visual magnitudes typical for globular 
clusters. To test fllS]) . we propose (1) resolving the kinetic energy in tidal tails outside the 
tidal sphere, (2) the velocity dispersion in the cluster and (3) comparing for the ratio of the 
two with the predicted value 

efc> 0.71 (l- 0.097y + 0(1/2)) Efc, (50) 

were equality is approached at large distances away from the cluster, when the gravitational 
binding energy to the cluster can be neglected. Following §2, R = l.SOrh in NE^ 
gives a mean (three-dimensional) velocity dispersion in the cluster. 



a 



0.19GM. 



(51) 
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and a scale for bipolar outflow 



Vt> 



0.84(l- 0.0492/+(y2))a 



(52) 



by ( l50l) . The velocity of the tidal fr ont may be l arger than V^, as it is determined by the 
fastest escaping stars. For instance, iLeon et al.l (120001 ) report for NGC 6254 a projected 
expansion velocity 7 km s~^ of tidal material at a distance of r =150 pc {D =4.1 kpc) and a 
projected velocity dispersion 6.6 km s~^ in the cluster, were it is noted that the gravitational 
binding energy to the cluster is negligible at this distance. Evidently, an expansion velocity 
exceeding the velocity dispersion in the cluster satisfies (!52l) . The correlation ( l50l) is relatively 
insensitive to various cluster parameters, such as the total mass, which otherwise is not 
well defined in view of an u ncertain but possible substantial population of white dwarfs 
teaumgardt fc Makinolboosh . 



In Table 1 we give an illustrate sample of our model predictions, where we assume a 
mean stella mass ffi = ^Mq to derive N = Mc/fh for each cluster. We focus on tidal tails 
created over the period T = TylO^ yr of a few times ten million years since the last disk 
crossing, to circumvent the uncertain fate of extended tidal tails shown in Fig. [2] following 
the most recent disk crossing. Scaled to one-quarter of the period P = Pq Myr of the orbit 
around the galactic center, the associated outflow of evaporating stars satisfies 



n 



PNf, 



N 



it 



0.275 X N^Pq 



relax 



rh 



109 yr 



(53) 



with N = 10^ N^. 



By angular momentum conservation, bipolar outflows emanating from the two Lagrange 
points Li ar id deflect into leading and trailing streams inside and outside the orbit of 
the cluster (iHeggid l200ll ) . We express it by a correlation of the radial separation to the 
temperature of the cluster according to an angular separation 



Ar 



a 



D 



3.4 



a 



0.55 ai^y^ deg. 



(54) 



where O = y ^ with Mg = 3 x lO^^M© ( Odenkirchen et al. 1997 ) denotes the orbital period 
of the cluster about the galactic center at a distance Dgc = -D4 10 kpc and cr = o"i 1 km s~^. 
The projected angular separation follows from the known positions of the globular clusters, 
shown in the last column of Table 1 in the approximation of circular motions. 
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9. Conclusions and observational outlook 



Stellar clusters are open systems that radiate stars and energy by evaporation, mostly 
by diffusion in momentum space subject to a singular perturbation by rare but finite jumps 
due to hard 2- or possibly 3-body interactions. We identify a parameter range for relaxation- 
limited evaporation applicable to clusters in their pre-coUapse phase (generally of KM mor- 
phology) with N » Nc and small to moderate tidal ratios y < 1, and develop a model for 
the evaporative lifetime for single-mass clusters. For a Maxwell velocity distribution that 
extends beyond the escape velocity, evaporation is bright as characterized by an inequality 
between the Kelvin- Helmholtz (ITS]) and Ambartsumian-Spitzer time scales satisfying 



txH < tAS 



(55) 



due to Jh > Jn, wherein evaporation is largely driven by energy-loss rather than mass loss. 
The evolution isolated clusters, as a reference to more realistic modeling, is hereby described 
by a finite-time singularity solution 



{N,R,T} oc (1 



r \ {0.0905, 0.6351, -0.5405} 
T ev^ 



(56) 



It shows that the temperature and the radius evolutions (l56l) are relatively less shallow than 
( !33l) . but still evolve only moderately away from the time of complete evaporation according 
to (l56l). Here, the indices 0.0905 and 0.6351 in (1561) are close to the values 0.085 and 0.695, 



respectively, observed in N-body simulations ( jHeggie fc HutI 120031 ). We do not attribute the 
relatively small t^v (rel ative to t^c) to a radial stratification of the relaxation time with 
/jy = as proposed in (I Johnstone! 119931). It also at odds w i th an older claim that jn < 
(with commensurably txH > ^as) by von Hoerner ( 1958 ): King ( 1958l lbl). Similarly, stars 
escaping by ejections from last scattering events from within the cluster should be contrasted 
with low velocity evaporation arising from diffusive escape through a low temperature outer 
boundary, commonly used in the Fok ker-Planck app roximations as if the cluster were opaque 
throughout up to its tidal boundary ( ISpitzerlll987l ). 



The agreement of ([2]) with numerical simulations supports the idea that energy loss in 
accord with ( 155|) is important. For a Maxwell velocity distribution, ([2]) is shorter than that 
based on particle loss alone by a factor of about three due energy loss 1^^. It would be of 
interest to i dentify the various s c attering process e s that may give r is e to a relative l y larg e 



fu. see e g. IPetrovskaval (Il970| ): iR.ettereii (Il979[ ): I.Tohnstond (Il993[ ): iMevlan et al.l (120011); 



Ashurovl (120041 ). and determine the implied shortening of the evaporative lifetime following 
§5.1. These studies fall outside the scope of the present work, however. Instead, we focus on 
some observational tests in 58. 
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The presence of a tidal field modifies the net rate of evaporation by two effects. It 
generally enhances evaporation by lowering the threshold on the kinetic energy for stars to 
spill over the tidal radius, but the escape probably is attenuated by the generally anisotropic 
tidal field as defined by the Lagrange points Li and L2 described by the grey body factor fl5^ . 
For moderate strengths of the tidal field, the resulting evaporation rate can be expanded in 
a regular perturbation in the tidal ratio ([1]). 

In studying the globular clusters of the Milky Way, we note that they (1) amply satisfy 
the criterion for relaxation-limited evaporation f ll6p . (2) their tidal rations satisfy y « 1, (3) 
most of them are of KM morphology indicative of a pre-core collapse state. For this reason, 
the perturbative expansion ([2]) for their remaining relaxation-limited evaporative lifetimes is 
believed to be applicable, and it is found to be consistent with their current ages. 

The energy of the escaping stars is directly correlated to the temperature of the cluster 
(!54|) . and the outflows thus produced develop a bi-polar 5"— shaped morphology in the plane 
of the Lagrange points Li and L2 and the orbit of the cluster. For those tidal tails where 
the S"— shaped tail can be resolved observationally, the de-projected orbital displacement of 
the tails is representative for the outflow velocity by the correlations 



If periodically destroyed by disk crossings, these tails hereby re-appear on a time scale 
of 10 Myr with a length scale of a few times the tidal radius for clusters that are not 
too small in mass and velocity dispersions. Detailed observational studies are proposed to 
test (1521) on a wide field of view around globular clusters using upcoming high resolution 
photometric and spectroscopic stellar survey instruments. The results may also be used for 
comparison with alternative mecha nisms for ind ucing tidal tails, such b y tidal heating or 



shocks in response to disk crossings (jHeggiell200ll ) and velocity anisotropy (ITakahashi fc Lee 



2OOOI ). We remark that destruction by tidal fields of an initial distribution might have been 
more severe around spiral galaxies than around the relatively larger elliptical galaxies, which 
tend to have relatively large numbers of globular cluster s in proportion to their luminosity 



(IHarris &: van den Berg Il98ll : iBurkert &: Tremaind 120101 ) . 



The upcoming large field of view photometric and spectroscopic instruments LSST, 
BigBOSS, GAIA and JASMINE seem to be ideally suited for detailed measurements on tidal 
tails in energy ( !52|) . stellar count (l53ll and morphology (15^ to advance an observational test 
of the present thermodynamic model. 
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TABLE I. Data on selected globular clusters derived from the catalogue of iHarrisI (jl996[ ) 
(revision 2010) and model predictions for the outflow velocity Vt from a given by (15T!) . the 
number of escaping stars n per one-quarter orbital period by ( !53l) and the projected angular 
separation ao between the leading and trailing tails inside and outside the orbit derived 
from using positional data. Where available, a comparison with observational data is 
included. 



Name" 








Vh 


Uh 


A 




a 


Vt 


n 






kpc 




pc 




Gyr 




pc 


km/s 


km/s 




arcmin 


NGC 1041 


4.5(7.4) 


1.45 


64 


0.065 


3.5 


3.7 


4.1 


16 


14 


140 


120(~ 100)^ 


Pal 41 


109(111) 


0.0541 


100 


0.162 


2.6 


2.0 


16 


1.6 


1.3 


161 


1 


Pal 51 


23.2(18.6) 


0.0284 


101 


0.182 


6.6 


2.1 


18 


1.1 


0.9 


6 


0.4(< 30)-^''^ 


Pal 14^ 


73.0(71.6) 


0.0200 


118 


0.219 


10 


1.5 


26 


0.8 


0.6 


23 


0.06 


NGC 51391 


5.2(6.4) 


2.64 


104 


0.073 


12 


1.2 


7.6 


16 


14 


64 


75(~ 100)-^ 


NGC 59041 


7.5(6.2) 


0.834 


52 


0.074 


2.5 


5.3 


3.9 


13 


11 


87 


43< 60)'' 


NGC 4590^ 


10.3(10.2) 


0.306 


89 


0.051 


1.9 


6.4 


4.5 


7.3 


6.1 


82 


22(< 60)^^ 


NGC 5466^ 


16.0(16.3) 


0.133 


162 


0.066 


5.7 


2.5 


11 


3.1 


2.6 


29 


4(< 60)3 


NGC 6254^ 


4.4(4.6) 


0.225 


30 


0.084 


0.8 


18 


2.5 


8.4 


7.0 


51 


70(~ 100)-^ 


NGC 6121^ 


1.7'^(5.9) 


0.225 


19 


0.112 


0.9 


23 


2.1 


9.0 


7.6 


93 


324 


NGC 67522 


4.0(5.2) 


0.364 


52 


0.042 


0.7 


17 


2.2 


11 


9.5 


96 


46 


NGC 7078^ 


10.4(10.4^ 


0.984 


62 


0.049 


2.1 


7.1 


3.0 


16 


13 


301 


8.7 



1 King morphology (jChernoff fc piorgovksilll989l) 

2 Post core-collapse morphology ( jChernoff fc Djorgovksilll989h 
NGC 104=47 Tuc, NGC 5139=w Centauri, NGC 5904=M5, 

6121=M4, NGC 6254=M10, NGC 7078=M15 



NGC 4590=M68, NGC 



Distances to the Sun (galacti c center) with uncertainties of 6% (IChaboyeii 120071 ) 



Chernoff fc Diorgovksil (119891 1 



Richer et al 



2004h 



^ Model predictions for projected angular separation between leading and trailing tidal 
streams (observed values, if known) from the position of the globular cluster assuming 
circular motion and M„ = 3 x IO^^Mq. 



Leon et al. (2000) 



Jordi fc Grebell (120101 ) 



Odenkirchen et al.l (120041 ): iBelokurov et all (120061 ): iFellhauer et al.l (120071 ) 



